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Steinhausen's hypothesis, that a simple torsion-pendulum model adequately describes the dynamics of 
semicircular-canal action, has received its sharp of criticism- This is primarily because of the fact 
that the time courses of* afferent discharges and vestibular reactions to angular-acceleration stimuli 
only approximately relate to the time course obtained as a solution to the trrsion-pendulum model for the 
same angular accelerations. Thus, there is a need t o determine a generalized transfer function whi^h 
delineates those components that provide meaning ti, Uteinhauspn's hypothesis from those which describe 
other phenomena such as "adaptation". 


INTRODUCTION 


Rotational head movements stimulate the sensory end organs of the semicircular finals of the inner ear, 
causing compensatory eye movements (nystagmus), appropriate to the plane of rotation, via a vestibtilo- 
ocular reflex arc. However, the sensory epithelia ..f the semicircular canals, the erlstae ampul lares, 
are relatively inaccessible for direct functional .study during normal physiological head movements. 

In *he first place, they lie deep within the skull, enclosed by spongy and petrous bone and the osseous 
labyrinth. Secondly, any technique which exposes the hair cells or their v>erv«* supply fo r neurotiological 
experimentation necessarily disrupts the vascular supply to the neuroepithelium and interferes with the 
ionic balance between the perilymph and endolymph. Thirdly, even when the membranous ampullae are 
exposed, it is difficult to perceive the extent to which the normally invisible cupula covers the 
crista and fills the ampulla of the duct. (It is the cupula which couples the mechanical energy of head 
motion to the receptor hair cells of the crista for transduction into neural impulses at- the afferent 
terminals.) Thus, it is not surprising that, for over a century, the dynamic response characteristics 
of the semicircular canals have not been completely elucidated despite continuous investigation by 
many vestibular physiologists. 

Historically, until the time of Fluorens (lA), the semicircular canals were associated with the 
perception of sound. By sectioning the membranous semicireti 1ar canals in the pigeon and the rabbit, 

Fluorens established that eye, head, and body movements always occur in the same plane as that of the 
injured canal, floltz (l6), also investigating the pigeon, first associated these canals with bodily 
equilibrium. Breuer (5), Crum Brown (9), and Mach (.’6) suggested that the crist.ae ampullares 
respond to rotatory accelerations of the head. More.,ver, they also hypothesized that the motion of 
the endolymph within the canal was re-sponnible for eliciting the neural activity in the nerve endings 
(hydrodynamic theory). (Mach, however, later appears to have favoured a hydrostatic theory, i.e., 
that pressure is the adequate stimulus—see refs. 7 f *, Pfi. ) Evald (li), using an "hydraulic hammer", 
produced artificial endolymphatic movements in the canals of pigeons and demonstrated that the 
direction of head nystagmus was dependent upon the direction of endolymph flow. Evald's "law" was 
placed on a firmer footing when Lovenstein and Sand (?5) established that the action-potential 
(impulse) frequency in fibers innervating the cristae varied according to the direction of angular 
movement of the head. 

In order to describe the dynamics of the semicire-ijar-canal system, rteinhausen ( 7 6, ^7), and, 
subsequently, F^mond, Oroen and .Tonkees (10), proposed the so-called torsion-pendul um model. In 
this model, the Instantaneous angular deflection of the cvjpula, E(t), ennped by an instantaneous 
angular acceleration of the head, o(t), is governed by the linear differential equation 


0 d + n - £ - ( ^1 ♦ Atft ) = YOa(t), 

dt ? dt 


(1) 


where 0 is the effective moment of Inertia of the endolymph in the semicircular canal; H is the 
viscous damping momeni per unit angular velocity of endolymph relative to the wall >T the membranous 
canal; A is the elastic-restoring moment per unit angular displacement of the end-lymph relative to 
the membranous ampullary walls; and t is a constant of proportionality relating ratio volumetric 
displacements (between the cupula and the canal endolymph). The transfer function pq, (l) can 
be written as: 
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; ne way t,. 1 ■ the adequacy •?' this model is t.> r*»c the neural activity <f the primary afferent 
fibers innervating t he semicircular canals a.; they resjiond t.<> appropriate stimulation. This has been 
done f<-r a variety -*f specie;: ,f animals (2, 1, 12, 17, IP, * ! 3, 2u, 33, 3l<, ■!>, 3ft) and the model has been 
found to he deficient ., The present set of experiments was designed t' _ study the mechanoneural response 
ohuruoterist lea of primary afferent, semicircular-canal units in the pigeon, in order to determine 
whether r not a generalized transfer function e,.uid be obtained which would also describe similar 
neurodynamic;: in other species. (In testing to this model, the assumption is made that the afferent 
activity is proportional to the cupular displacement.) 


METHODS 

White King pigeons ( Oolumba H vi a) were surgically prepared for microelectrode recording from 
peripheral units which innervated the semicircular canals (primarily, the anterior semicircular canal). 
The animal (with head immobilized) was oriented on board a rotatory device so that the center of its 
head plane was coincident with t.he plane of rotation (according to the type of canal being recorded 
from). 

The main rotatory sequence consisted of a scries of sinusoidal angular accelerations, 

aft) = a sin Pirft, (3) 

m 

which were delivered to anesthetized preparations at frequencies, f, from 0.01 to 10 Hz with peak 
angular accelerations, a m , of 2.0, 1».0, 8.0, 12.0, and 20,0 degrees/s^ (t ^ real time in Eq. (3)). 

The single unit neural activity was amplified, displayed on an oscilloscope and recorded on magnetic 
tape according to conventional techniques. (One channel of the tape recorder was used for voice 
commentary; another for the stimulus reference signal.) 

The stimulus reference signal was used to trigger a physiological signal analyzer (Nicolet 
Instruments Inc.) t.o count and store the taped, entrained action potentials in preselected appropriate 
time periods (bins) for up to U096 sequential bins (see Fig. 1 for typical binned response). 

Fourier techniques were used on the binned neural data to obtain the magnitude of the peak amplitude of 
the response and the temporal (phase) relationship between the angular acceleratory stimulus and the 
fundamental component of the neural response. 

Amplitude- and phase-values were usea in a curve-fitting program on a PDF'-ll/bO minicomputer 
(Digital Equipment Corp.) to provide a distinct mathematical expression for the best-fitting transfer 
function for a linear system. (The steady-state frequency response Tor a linear system to a sinusoidal 
input angular acceleration may be found from gain, 

!0(f)| = ((He 0(f)) L> + (Im (Hf))' - ’} 5 *, (M 

and phase, 

0(f) = tan -1 dm <i(f)/Re G(r)), (5) 

spectra (Bode plots—see ref. 30), where Re 0(f) and Im 0(f) are the real and imaginary parts of G(f), 
respectively.) 

The method employs techniques in non-linear, least squares approx.mations and is applicable to 
both high- and low-order transfer functions 0»). Ori the final Interpol* tlon in the program, the least 
squares error, LSE, of the best-fitting transfer function is obtained as 

L&E * Y T Y, (<>) 

where Y is a residuals matrix (Y*^= transpose of Y) which is made up of error terms that express the 
differences between the experimental data and the model which is to be fitted. The mean square error, 

M15E * LSE/(2L - S), CO 

compares the goodness-of-fit of the derived transfer function to that of the experimental data (2L £ 
number of real and imaginary components of the data points, and i number of parameters in the transfer 
function). 

A more detailed description of the methodology may be found in a special monograph (7) anu elsewhere 
(8, 21). 


result: 4 and discussion 

The simplest transfer function that fitted the data for all units ia of the form. 
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( 0 ) 


where s i« Cract i1 — -rder differential operator with 0 < k < 1, and n is a gain constant with 
units in impulu*-:--s~ 1 /ii» grct-:-, (.’1). The Bode pints for four of the units are shown in Figs. 2-5, 

together with their Best-fitting ami the best-fitting torsion-pendulum model 


(In rigs. 2-5, the effects of the To — 2.0 ms ( vide infra ) contributes no more than ll to the gain spectra 
between f - 0.01 and 10 Hz.; consequently, the single-pole transfer function, 0"(s), was an adequate 
representation of the torsion-pendulum model.) As is evident i» the plots in Figs. 2-5, G'(s) is a 
much be* ter *'it to the data than is G"(s) (of. the M3F. values for the two models); and, in particular, 
the fit appears to improve with increasing k. 

K 

Wnat, then, is the significance of k, or better still, s ? Elsewhere, it is shown that 


M 

* K n 

i = l 


( 10 ) 


where K and Bj are constants, H s * (in theory, but finite when fitting Eq. (10) via a digital 
computer), and and tj j {Cq : 0) are time constants (21), Interestingly, when M = 1, Eq. (10) 
becomes 


( 11 ) 


which has previously been defined as the transfer function of the adaptation operator (27, 1*0). Thus, 
appears to be a form of adaptation. Work by Thorson arid Riederman-Thorson (39) suggests that 
s k represents a distributed relaxation process which is inherent in the sensory-adaptation mechanics 
of Li mu 1 us photoreceptors, vertebrate retinal receptors, chemoreceptors, and other mechanoreceptors. 
Investigations by Taglietti, Rossi and Casella (30) further suggest that s* 1 likely represents a 
relaxation phenomenon consisting of a time-varying intracellular electrogenic process, the com¬ 
ponents of which are summed with the generator potential in the receptor hair cell. 


The coefficient, of variation, CV, was determined as the ratio of the standard deviation of intervals 
to the mean interval, as obtained from int.erspike-interval distributions of spontaneous single unit 
activity. When a regression of CV on k was made for 20 units, a statistically-signifleant product- 
moment correlation (r = 0.301*, F’ < 0.05) was obtained (21). Thus, the larger the CV is, the larger the 
value of k and, consequently, the amount of adaptation. Other work by Goldberg and Fernandez (15). in 
squirrel monkeys, shows that the CV Is statistically correlated with semicircular-canal afferent 
fiber conduction rates. The thicker fibers have faster conduction rates and larger CVs. Together, 
these findings suggest that sensor 1 adaptation phenomena are directly dependent on the innervation 
pattern of the afferent fibers. 

The transfer function, G*(s), differs also from that of Ci"(s) in that is not single-valued as 
it is in the torsion-pendulum model; rather, it is unit dependent, taking on values from t^ * 1*.1*5 to 22.17 
s (mean ± GEM * 10.2l» ♦ 1.20 s) (21). (in fitting Eq. (0), the coupling between k and would account 
for some of the five-fold range of values that were determined for t L . However, there is sufficient 
indication from other studies (32) that the response dynamics of small groups of contiguous hair cells 
are quite different from those of other groups.) Realizing that the nair-cell tufts are stiff (lU), 
that their lengths vary according to their position on the sensory epithelium (crista) (22), that the 
number and thickness of the stereocilia can be variable (22), and that the mechanical properties of the 
cupula are not necessarily uniform across the crista (29), then it is plausible that could have 
a regional distribution. 


The form 

W 1 C« K 

(x^nHigS* 1 ) • (l2) 

or sometimes, 

G(s)(t s+1) (13) 

m 

was fitted to published afferent-response data in the squirrel monkey ( Salmiri sciureus ) (12), the 
frog (Rana esculents and R. temporarla )(3), the gerbil ( Meriones ungulculatusT (35), and the guit&rfish 
( Rhlnobatos productus ) (33). The parameter x m = high frequency time constant which results from both 
the displacement and the rate of displacement of the cupula. Such a term has been obtained from 
analysis of vestibular-driven eye movements In man (Benson and Sternfeld, cited in ref. l), and 
primary afferent canal responses in the squirrel monkey (12) and the elasmobranch fish (2U). Table 1 
lists k, T., t, and for these four species and for unit responses to white-noise stimuli obtained 
from pigeons that were primarily enc^phale isol£ preparations. For all species listed in Table 1, the 
MGE using Eqs. (12) or (13) was comparable to or smaller than that obtained with other models. Further 
details are given elsewhere (0). 




! .lift.i >ti (Eqs. M:'> nr (13)) describing jemirircular-ranaJ 
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TABLE 1: Parameters >. f a genera 
dynamics in five self 1 




' '• T m 
s P* cies i ms 1 (s) 

frequency range, 
f • H? ) 


Squirrel monney 

(a) "regular" ur.i’. 

(b) "irregular" «. •* 

Gerbi I 

(a) CV < 0 1 

(b) CV > 0. I 

Fmg 

Guitarfi sh 


0.0125-R 

0.0125-8 


0.01-5 

0.01-5 


0.0125-0.5 

0 ,w2-U 


0.5-16 


In Table 1, the frequency ■••» 
adequate value of Tg. In gene/ ■» 
from the biophysical proper* i«v 
of the membranous vestibular at t *» 
Money and colleagues (U) f r *;• 
(f = l/2tiTj,) of 80 Hz in the tv* 
from canal' afferents in the g . : 
given evidence that 5.70 ^ t • 
fitting to afferent data fn t k <«* 


* » restricted in the frog and guitarfish to utilize an 

• • (s 0/11) in the fits have been predetermined indirectly 

.ymph and the dimensions of the pertinent anatomical features 

• .*a.>Jt* irj = 2 ms, which was determined biophysically by 
. .cgf-it; that there should be an upper break t'requency 

r«'..pi ii.-.e dynamics in semicircular-canal afferents. Recording 
. .t. ,. n'k frequency up to 70 Hz, Hartmann and Klinke (l8) have 
•t id the same range as those determined empirically for 

i. , .’1, <5). 
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